Interaction-induced shift of the cyclotron resonance in graphene via infrared spectroscopy 
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We report a study of the cyclotron resonance (CR) transitions to and from the unusual n — Landau level 
(LL) in monolayer graphene. Unexpectedly, we find the CR transition energy exhibits large (up to 10%) and 
non-monotonic shifts as a function of the LL filling factor, with the energy being largest at half-filling of the 
n — level. The magnitude of these shifts, and their magnetic field dependence, suggests that an interaction- 
enhanced energy gap opens in the n — level at high magnetic fields. Such interaction effects normally have 
limited impact on the CR due to Kohn's theorem [W. Kohn, Phys. Rev. 123, 1242 (1961)], which does not apply 
in graphene as a consequence of the underlying linear band structure. 

PACS numbers: 78.66.Tr; 71.70.Di; 76.40.+b 



The intriguing electronic properties of graphene in a strong 
magnetic field were first demonstrated by the observation of 
an extraordinary 'half-integer' quantum Hall effect CHS]], re- 
sulting from the existence of an unusual n = Landau level 
(LL). Later evidence from transport measurements indicates 
that the four-fold degeneracy of this n = level is completely 
lifted in high magnetic fields OH, and a gapped state forms 
I8UTT1I. This state and the origin of degeneracy-breaking in 
the n = level is a matter of intense theoretical study, with 
electron-electron interactions expected to play a critical role 
fT2l . However, since the charge transport in earlier experi- 
mental works is either local, or dominated by the sample edge, 
direct access to the intrinsic properties of graphene in the bulk 
remains limited. Additionally, the diverging resistance near 
half-filling of the n = level prohibits investigation of gapped 
states by thermal activation experiments. In this Letter, we 
report a study of the n = LL utilizing infrared (IR) mag- 
netospectroscopy, which is sensitive to the cyclotron orbits 
of charge carriers throughout the entire graphene sheet. As 
the LL filling factor is changed, we observe unexpected and 
sizable shifts in the cyclotron resonance (CR) transition ener- 
gies. We interpret this as due to electron-electron interactions 
which create a gap in the n = LL, thereby affecting the 
energies of CR transitions to and from this level. 

In a perpendicular magnetic field, B, the LL structure of 
graphene has an unusual energy dependence on both B and 
the LL index, n, described by |[T3lfT5l 



(1) 



Here h is Planck's constant, 
locity, and Ib 



c ~ 1 x 10 6 m/s the band ve- 
the magnetic length. The index 



n runs over both positive and negative integers representing 
electron and hole-like states, and includes a field-independent 
n = LL which is located at zero energy. This level, in com- 
bination with a four-fold level degeneracy due to the two-fold 
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FIG. 1: (a) Optical microscope image of sample A. (b) Measured 
Hall (blue) and longitudinal (red) 4-wire resistances of sample A. 
(c) Calculated Hall (black) and longitudinal (green) conductivities, 
showing a Hall plateau forming at v — 0. (d) Schematic of ex- 
perimental setup showing IR light focused through graphene. The 
varying bolometer resistance tracks changes in the transmitted light. 



electronic spin and pseudo-spin (valley) degeneracies, leads 

to a unique sequence of LL filling factors, v = ±2, ±6 

Here v = 2i\n s l 2 B , where n s is the carrier density. In con- 
trast to conventional 2D systems, in which even the lowest LL 
experiences an energy shift with B, in graphene the n = 
LL remains put at zero energy, is shared by the conduction 
and valence bands and, at charge neutrality for vanishing is, is 
filled half with electrons and half with holes. 

Motivated by recent experiments trjlfTTlL we have per- 
formed a careful study of LL transitions to and from the n = 



level via CR [16-21']. The five samples, A-E, used in this work 
consist of monolayer graphene exfoliated from Kish graphite 
onto Si/Si02 wafer substrates (22), with Cr/Au contacts de- 
fined by standard thin film processing for in situ measurement 
of sample transport properties. Sample sizes and mobilities 
(at a hole density of 1.7xl0 12 cm -2 ) range from 330-1,100 
/mi 2 and 3,000-17,000 cm 2 /(Vs), respectively. The lightly- 
doped Si substrates transmit mid-IR light, but remain suffi- 
ciently conductive for use as a back gate to vary the carrier 
density. Samples are mounted at the focus of a parabolic mir- 
ror and cooled via He exchange gas in a probe inserted into 
a liquid He cryostat. Broadband IR light from a Bruker IFS 
66 Fourier transform spectrometer reaches the graphene via 
direct optics, and the transmitted light is detected by a com- 
posite Si bolometer placed immediately below the sample; see 
Fig. 1 (d). The samples are wedged to suppress Fabry-Perot 
interference. We determine the CR spectrum at constant B 
field as in Ref. ifTTl . recording IR transmission through the 
sample at two gate voltages, V g and Vt,. The change in trans- 
mission, AT/T = 1 — T(V)/T(V b ), is measured over a range 
of V g (y). The background scans, T(V b ), are at v = ±10. 

In Fig. 1 (a) we show an optical microscope image of our 
highest mobility monolayer graphene device, sample A, with 
mobility /i = 17, 000 cm 2 /(Vs). The Hall and longitudinal 
resistances for this sample, and the calculated conductivities, 
are shown as a function of V g in Fig. 1 (b) and (c) respectively. 
A plateau in the Hall conductivity is clearly evident near v — 
0, as in Ref. 0013. 

In Fig. 2 (a) we plot the change in IR transmission, AT/T, 
through sample A for a range of v at a field of B = 18 T. From 
bottom to top, the filling factor decreases from v = -f 4 x /2 to 
v = —3^2 in steps of Au = — V 2 - Physically, this amounts 
to shifting the Fermi level, Ep, through a fixed LL spectrum 
starting just above half-filling of the n = 1 LL and ending 
below half-filling of the n — — 1 LL. The broad peaks corre- 
spond to the resonant excitation of carriers between the lowest 
LLs, via the n = — 1 — » and/or the — » 1 transition. The 
sharper spikes at 157 meV and 196 meV are due to harmon- 
ics of 60 Hz electrical noise. Additional small features resem- 
bling peaks or shoulders on the lineshape are visible in several 
traces. These are suggestive of splittings of the CR line. Such 
an observation would be intriguing and underscores the need 
for higher quality samples, but here we focus on the main CR 
peak energies traced by the dashed line in Fig. 2 (a). 

Closer inspection of the CR traces reveals a distinct asym- 
metry in the lineshape with the high-energy side falling off 
more slowly than the low-energy side. This distortion of 
the Lorentzian resonance originates in multiple reflections be- 
tween the graphene and the wedged Si/SiC>2 substrate l23l . 
To extract the physical parameters for each resonance, we 
have performed detailed curve fits by modeling the transmis- 
sion through a multilayer system including substrate phonon 
modes and multiple reflections/refractions at the graphene- 
SiC>2 and SiCVSi interfaces. The dielectric properties of the 
substrate used in the fits are determined experimentally by a 
combination of IR spectroscopy and ellipsometry l24l . Cal- 
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FIG. 2: (a) The change in transmission, AT/T, for CR traces in 
monolayer graphene at fixed field B = 18 T over a range of v. 
Traces are offset for clarity. The bottom trace is v — 4 1 /2\ v of 
successive traces decreases by 1 /2. The dashed line is a guide to the 
eye, marking the peak energies determined by curve fits, (b) Exam- 
ple of a fit (solid red line) to the CR line shape at v — I 1 / 2 using a 
multilayer IR transmission model, (c) CR peak energies plotted vs. 
v. The area of the markers represents the relative intensity. 



culated fits are compared to the data in successive iterations, 
as the model resonance energy, intensity and width are varied 
until best fits are achieved as determined by inspection. An 
example is shown in Fig. 2 (b) for the v — +1 1 / 2 trace. We 
note the peak energies from this fitting fall at energies slightly 
below the intensity maxima l25l . 

Far from remaining constant as the occupancy of the n = 
LL is changed, the CR exhibits a marked non-monotonic 
dependence on v, resulting in a nearly symmetric "W"-shape 
as shown in Fig. 2 (c). Two sharp minima in the CR energy 
occur at v = ±2, with a more rounded maxima at v = that 
is nearly 10% higher than the values at \v\ = 2. We note the 
increase in energy occurring for \v\ > 2 is accompanied by a 
loss in intensity, with the peaks fading away beyond \v\ =4. 
This is due to the onset of the next LL transition, n = 1 — >• 2 
or —2 — >> — 1, as E F is shifted farther from charge neutrality 
and either the n = +1 or — 1 LL becomes partially occupied. 

The considerable shifts seen in Fig. 2 (c) are not reflected in 
the interband transitions that are excited simultaneously, albeit 
at higher energies. In Fig. 3 we show the CR data for v = 
±2 and in sample B over an extended energy range which 
also includes the first degenerate pair of interband transitions, 
n = — 2 — )> 1 and — 1 — »• 2. The transitions responsible for 
the two resonances in each trace are illustrated schematically 
in the upper left inset to Fig. 3. The interband peak (at 440 
meV) shifts by less than 1% as v is changed. In contrast, and 
consistent with the behavior seen in Fig. 2, the lower energy 
peaks (at 165-180 meV) show the same behavior as in sample 
A, with a strong upshift (« 8%) in energy at v = 0. 

The increase in CR energy of the lowest LL transitions at 
v = depends on both B and the sample quality. Figure 
4 shows the B field dependence in all five samples at fields 
up to 3 1 T. Here we characterize the overall magnitude of the 
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FIG. 3: Two CR peaks visible in AT/T traces from sample B, at 
B = 18 T. The underlying transitions are n = — ► 1 and/or n = 
— 1 — )► near 170 meV (blue arrows), and n — — 1 — >> +2 and 
n — —2 — )■ +lat 440 meV (red arrows). Traces are offset for 
clarity. Left inset: LL schematic showing transitions at v — +2, 
when Ef lies between the zeroth and first LLs. Right inset: proposed 
LL schematic at charge neutrality {y — 0), showing a gap in the 
n — LL and the corresponding lowest energy CR transitions. 



shift, AE, as the difference in CR energies at v = +2 and 
v = 0. Open symbols show data from the highest mobility 
sample, A, and filled symbols are data from samples B-E. A 
clear trend of increasing AE (up to 23 meV) with increasing 
B emerges. We find AE is correlated with device quality: 
sample A shows the largest shifts, while lower mobility sam- 
ples have consistently smaller AE. Although the broader res- 
onances in the lower mobility samples lead to wide error bars, 
in sample A the uncertainty is sufficiently small to explore a 
functional dependence between AE and B. In the inset to 
Fig. 4, we plot the CR energies of the lowest LL transitions 
vs. y/B at two fields for sample A. As suggested by Eq. 1, 
these data are well described by straight lines constrained to 
pass through zero and differing only in slope, implying that 
a \fB dependence is followed in both cases. Further experi- 
ments in higher quality samples are needed to determine the 
precise field dependence of AE. 

We now discuss possible origins of the observed CR shift 
at v = 0. In the basic model of non-interacting carriers, the 
CR energy is the difference of the initial and final LL energies 
from Eq. 1. For \v\ < 2, two degenerate transitions occur, 
n = and 1. For 2 < v < 6 then = 1 2 also 

occurs, but with an energy reduced by a factor of y/2 + 1 that 
lies outside our experimental window. Thus we expect to see 
only a single CR peak with energy independent of v, which 
clearly is not observed. 

Disorder certainly plays a role, as seen in Fig. 4. We note a 
simple model where LLs are broadened only by short-ranged 
scattering predicts an upshift in the CR for half-filled LLs 
1261 , similar to our observation. However, one would expect 
such a CR shift to increase in more disordered samples, in 
contrast to the observed increase in AE with increasing mo- 
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FIG. 4: B field dependence of the CR shift, AE, between v — +2 
and 0, as a function of B for five samples. Open symbols are from the 
high mobility (/x = 17, 000 cm 2 /(Vs)) sample A resonances shown 
in Fig. 2; filled symbols show lower mobility samples (/i < 10, 000 
cm 2 /(Vs)). The dotted line is the bare Zeeman splitting, 2\ibB. Er- 
ror bars are estimated by varying the best fit parameters. Inset: v — 
and +2 CR energies of sample A at two fields, plotted against y/B. 
Dashed lines are least-squares linear fits through zero. 



bility. Also, the disorder potential in our samples is likely 
dominated by charged impurities (27) . Screening may affect 
the CR transitions as well. However, we expect increased 
screening to soften the CR |28], making for an M-shaped v- 
dependence rather than the W-shape of Fig. 2 (c); or else to 
narrow the resonance linewidth in half-filled levels [ 29 ] which 
runs counter to the bulk of our data. In this context, we note 
the broadening may be due to the linear dispersion alone (30). 

Single particle spin-related effects are also an unlikely ori- 
gin for our observations. Even at 3 1 T the bare Zeeman energy 
amounts to only 3.6 meV, several times smaller than the ob- 
served shifts (see Fig. 4). Also, as spin is conserved in the 
CR transition, simple Zeeman effects should not be observed. 
There remains the possibility of a ^-dependent g-factor en- 
hancement. Large enhancements have been seen in other 2D 
systems (3D . However, these appear only in transport data 
with the interpretation invoking many-body effects. 

The linear dispersion of graphene implies the CR can be in- 
fluenced by many body effects (32). In our experiments, the 
clear dependence of the CR energies on v- and hence E^- 
and the magnitude of the shifts suggest that interactions at 
the Fermi surface are an important ingredient. The contribu- 
tion of electron-electron interactions to the dispersion of low- 
lying inter-LL excitations has been calculated for graphene in 
Ref. [33], and also in Ref. (32) where the focus is restricted 
to those modes that contribute to CR. In these works, inter- 
actions are generally found to increase the excitation ener- 
gies. In particular, the calculations show that when v > 2 
(y < -2), the n = — )> 1 (n = -1 0) transition de- 
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velops a splitting where the higher energy peak continues to 
increase in energy with increasing (decreasing) v. Qualita- 
tively, this can describe the CR shifts we observe at \v\ > 2, 
if we assume the splitting is masked in our data by the broad- 
ened lineshape. On the other hand, neither work predicts a 
change in the energy of optically active modes for \v\ < 2, 
and so fail to capture the large upward shift observed in our 
CR data near v — 0. Yet in spite of this, interactions seem to 
be the origin for our observations: although at 18 T the bare 
Coulomb energy e 2 /(els) ~ 90 meV is much larger than the 
shifts we observe in Fig. 2, no other energy scale is available. 

Zhang et al have observed quantum Hall plateaus at v = 
and ±1 which they interpret as interaction-induced splittings 
of the n = level |6|. Naively, one could view our results as 
the consequence of such a splitting, because when the Fermi 
level passes through the n = LL, the presence of a gap leads 
to shifts in the transition energies to neighboring LLs; see the 
right inset to Fig. 3. This picture, while certainly an oversim- 
plification, is nonetheless appealing: it easily accounts for the 
v- and independence of the CR energies, and is consistent 
with the fact that only those transitions involving the n = 
LL are strongly affected. 

As to similar experiments, the CR data on epitaxial 
graphene |2T1 have not yielded similar behavior. Previous 
experiments on GaAs and Si/Si02 2D systems have shown 
unexpected shifts and splittings of the CR bearing some re- 
semblance to our data |[34H36l . However in detail the com- 
parisons break down. An increased effective mass of the car- 
riers (decreased CR energy) seen for full LLs in 2D GaAs 
was attributed to a particular impurity doping, to band non- 
parabolicity, or to ^-dependent screening, all of which are un- 
likely to be the source of our observations in graphene. The 
origin of the anomalous CR peaks and splittings in the Si data 
l36l was never resolved, but in general these features had no 
B dependence, with higher quality samples showing weaker 
anomalies, both in contrast to our data. 

In conclusion, we have observed sizable shifts of the CR 
energy as a function of the LL filling factor and applied B 
field for transitions involving the n — level in monolayer 
graphene. We suggest our results indicate the opening of an 
interaction-induced gap in the bulk of the graphene at v = 
which increases with increasing B field strength, and is de- 
tectable via CR as Kohn's theorem does not apply in graphene. 
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